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We report on a search for 5° d(*)o j^(*)o decays based on 85 x 10® BB events collected with 
the Belle detector at KEKB. The D"K^^ and B" D"K*^^ decays have been observed for 

the first time with the branching fractions B{B° D°K°) = (5.0li:2 ± 0.6) x lO"'' and B(B" 
D°K*^^) = (4.8ti;J ± 0-5) x 10~^ No significant signal has been found for the B° D*°K^*^° and 
B" ]j(*)Oj^*o (jecay modes, and upper limits at 90% CL are presented. 
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Since the recent discovery of CP violation in the B me- 
son systenijthrough the measurement of non-zero values 
for sin 2(j)i [l| , attention has turned towards the measure- 
ment of the other Unitary Triangle angles. Such measure- 
ments will allow tests of the Kobayashi-Maskawa ansatz 
and of the Standard Model. Precise measurements of the 
branching fraction for B° D^K*°, B° D°K*° and 
B° D°(jpK*° decays, where D^jp denotes D° or de- 
cay to a CP eigenstate, will allow a measurement of the 
angle 4>3 01]' The decay B^ D^K^ can also be used to 
measure time-dependent CP asymmetry in B decays 0] . 
So far no experimental information is available for any of 
these decays. 

In this Letter we report on a search for the B^ 
D(*)OK°, B° ^ D(*'>°K*° and B° ^ d(*'>°K*° de- 
cays with the Belle detector at the KEKB asymmet- 
ric energy e+e" collider [^. The results are based on a 
78 fb^^ data sample collected at the center-of-mass (CM) 
energy of the T{AS) resonance, which contains 85 x 10^ 
produced BB pairs. 

The Belle detector has been described elsewhere 0. 
Charged tracks are selected with a set of requirements 
based on the average hit residual and impact parameter 
relative to the interaction point (IP). We also require 
that the transverse momentum of the tracks be greater 
than 0.1 GeV/c in order to reduce the low momentum 
combinatorial background. 

For charged particle identification (PID), the com- 
bined information from specific ionization in the central 
drift chamber (dE/dx), time-of- flight scintillation coun- 
ters (TOP) and aerogel Cerenkov counters (ACC) is used. 
At large momenta (> 2.5 GeV/c) only the ACC and 
dE/dx are used. Charged kaons are selected with PID 
criteria that have an efficiency of 88%, a pion misidenti- 
fication probability of 8%, and negligible contamination 
from protons. All charged tracks having PID consistent 
with the pion hypothesis that are not identified as elec- 
trons are considered as pion candidates. 

Neutral kaons are reconstructed via the decay 
TT+TT" with no PID requirements for these pions. 
The two-pion invariant mass is required to be within 
6 MeV/c^ (~ 2.5a) of the nominal K'^ mass and the dis- 
placement of the TT'^TT~ vertex from the IP in the trans- 
verse (r — (p) plane is required to be between 0.2 cm and 
20 cm. The direction from the IP to the tt+tt" vertex is 
required to agree within 0.2 radians in the r — (f) plane 
with the combined momentum of the two pions. A pair of 
calorimeter showers not associated with charged tracks, 
with an invariant mass within 15 MeV/c^ of the nomi- 
nal 7r° mass is considered as a 7r° candidate. An energy 
deposition of at least 30 MeV and a photon-like shape 
are required for each shower. K*'^ candidates are recon- 
structed from K~'K~^ pairs with an invariant mass within 
50 MeV/c^ of the nominal mass. We reconstruct D° 
mesons in the decay channels: if~7r+, i4r~7r+7r~7r+ and 
K~'K~^'K^, using a requirement that the invariant mass 



be within 20 MeV/c^, 15 MeV/c^ and 25 MeV/c^ of 
the nominal D'^ mass, respectively. In each channel we 
further define a D° mass sideband region, with width 
twice that of the signal region. For the tt" from the 
D'^ i^~7r"'"7r° decay, we require that its momentum 
in the CM frame be greater than 0.4 GeV/c in order to 
reduce combinatorial background. mesons are re- 

constructed in the Z?*° D^tt^ decay mode. The mass 
difference between and I?" candidates is required to 
be within 4 MeV/c^ of the expected value (~ 4cr). 

We combine candidates with or K*° to form 

B mesons. Candidate events are identified by their CM 
energy difference, AE = (X^j^i) ^ ^b, and the beam 
constrained mass, Mbc = \/E'^ — {^iPiY, where E'b is 
the beam energy and pi and Ei are the momenta and 
energies of the B meson decay products in the CM frame. 
We select events with Mbc > 5.2 GeV/c^ and |A£;| < 
0.2 GeV, and define a B signal region of 5.272 GeV/c^ < 
Mbc < 5.288 GeV/c2 and |A£:| < 0.03 GeV. In the rare 
cases where there is more than one candidate in an event, 
the candidate with the and ^(*)° masses closest 

to their nominal values is chosen. We use Monte Carlo 
(MC) simulation to model the response of the detector 
and determine the efficiency 0. 

To suppress the large combinatorial background domi- 
nated by the two-jet-like e+e" qq continuum process, 
variables that characterize the event topology are used. 
We require |cos6'thr| < 0.80, where 0thr is the angle be- 
tween the thrust axis of the B candidate and that of the 
rest of the event. This requirement eliminates 77% of 
the continuum background and retains 78% of the sig- 
nal events. We also construct a Fisher discriminant, T, 
which is based on the production angle of the B can- 
didate, the angle of the B candidate thrust axis with 
respect to the beam axis, and nine parameters that char- 
acterize the momentum flow in the event relative to the 
B candidate thrust axis in the CM frame |^ . We impose 
a requirement on that rejects 67% of the remaining 
continuum background and retains 83% of the signal. 

Among other B decays, the most serious background 
comes from B° D-tt+, -> K^*'>°K-, K^*^°K-Tr°, 
K(*}°K-TT-7r+ and ^ D- K+ , D- K'^*'>°Tr~ , 
A'(*)°7r^7r°, ^(*'"7r^7r^7r+. These decays produce the 
same final state as the B'^ £)(*)o^(*)o gjgj^ai^ ^nd 
their product branching fractions are up to ten times 
higher than those expected for the signal. To suppress 
this type of background, we exclude candidates if the 
invariant mass of the combinations listed above is consis- 
tent with the D~ hypothesis within 25 MeV/c^ (~ 3a). 
The B° D*+K-, D*+ D°tt+ decay can also pro- 
duce the same final state as the 5° D°K*° decay. But 
this decay is kinematically separated from the signal; the 
invariant mass selection criteria for K*'^ candidates com- 
pletely eliminates this background. Another potential 
BB background comes from the B^ — > D^*^^p^ decay 
channel j^] with one pion from p'^ decay misidentified as 
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FIG. 1: AE (left) and A4"bc (right) distributions for the — > /)0^{*)o candidates. Points with errors represent the experimental 
data, hatched histograms show the Z)" mass sidebands and curves are the results of the fits. 



a kaon. The reconstructed K~tt~^ invariant mass spec- 
tra for these events overlap with the signal K*'^ mass 
region, while their AE distribution is shifted by about 
70 MeV/c^. We study this background using MC simu- 
lation. The contribution to the B° ~* D^*^°K*° signal 
region is found to be less than 0.2 events. We examined 
the possibility that other B meson decay modes might 
produce backgrounds that peak in the signal region by 
studying a MC sample of generic BB events that corre- 
sponds to about 1.5 times the data statistics. No other 
peaking backgrounds were found. 

The AE and A/bc distributions for B° -» D°K^*'>° can- 
didates are presented in Fig.m where all three decay 
modes are combined. Also shown in Fig. ^ by hatched 
histograms are the distributions for events in mass 
sideband. The sideband shape replicates the background 
shape well, confirming that the background is mainly 
combinatorial in nature. Clear signals are observed for 
the D°K° and D°K*° final states. As an additional cross 
check, we also study the Kg candidates' invariant mass 
and flight distance distributions and K*'^ candidates' in- 
variant mass and helicity distributions for these decays. 
The helicity angle 9k* is defined as the angle between 
the K**^ momentum in the B meson rest frame and the 
K~ momentum in the K*^ rest frame. The distributions 
mentioned above are shown in Fig. [21 where points with 
error bars are the results of fits to the AE spectra for 
experimental events in the corresponding bin, and his- 
tograms are signal MC. All distributions are consistent 
with the MC expectation. 

For each decay mode, the AE distribution is fit- 
ted with a Gaussian for signal and a linear function for 



background. The Gaussian mean value and width are 
fixed to the values from MC simulation of the signal 
events. The region AE' < —0.1 GeV is excluded from 
the fit to avoid contributions from other B decays, such 
as i? ^ _D(*)°^(*)°(7r) where (tt) denotes a possible addi- 
tional pion. For the Mbc distribution fit we use the sum of 
a signal Gaussian and an enipirical background function 
with a kinematic threshold [ini |. with a parameter fixed 
from the analysis of the off-resonance data. For the cal- 
culation of branching fractions, we use the signal yields 
determined from the fit to the AE distribution. This 
minimizes a possible bias from other B meson decays, 
which tend to peak in Mbc but not in AE. The fit re- 
sults are presented in TableHl where the listed efficiencies 
include intermediate branching fractions. The statistical 
significance of the signal quoted in Table is defined as 
•\/— 2 lii{£.o / Cmax) , where Cmax and Cq denote the max- 
imum likelihood with the nominal signal yield and the 
signal yield fixed at zero, respectively. 

For the final result we use a simultaneous fit to the AE 
distributions for the three 13° decay channels taking into 
account the corresponding detection efficiencies. The 
normalization of the background in each D'^ sub-mode is 
allowed to float while the signal yields are required to sat- 
isfy the constraint Ni = Ngs ■ B{B° -> D^*'>° K^*'>°) ■ Si , 
where the branching fraction B{B° Z)(*)0/i:(*)oj jg ^ 
fit parameter; Ngg is the number of BB pairs and is 
the efficiency, which includes all intermediate branching 
fractions. 

The statistical significances for the B'^ D^K'^ and 
B^ — > D'^K*'^ signals are higher than Scr. The sig- 
nals in the D*^K^*'^^ channels are not significant 
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FIG. 2: From left to right: K% candidates' invariant mass and flight distance for the D^K^ channel, isT*" candidates' 

invariant mass and helicity distributions for the — > D"K*'^ channel. 



TABLE I: Fit results, efficiencies, branching fractions and statistical significances for — > D^*'^ K^*''^ decays. 



Mode 






AE yield 


Mbc yield 


Eflicicncy (10"'^) 


B (10-5) 


Significance 


BO ^ 




* K--K+ 


9 3+4.2 


7 1+4.0 
'•^-3.4 


2.50 


4.4l2;0±0.5 


3.1a 


SO ^ 


DOrO, dO - 


> K~ H'^ Tl" Tl^ 


14.9l|;5 




2.36 


7.4l|;^ ±0.8 


3.6(T 


SO ^ 


DORO, dO - 


* K-n+nO 


8.7tt-j 




2.52 


4.0l2;2 ± 0-4 


1.9cr 


SO ^ 


DOrO^ simultaneous fit 


3l.5l?;2 


27 0+^-'' 


7.38 


5.0li;2 ±0.6 


5.1cr 


BO ^ 
BO ^ 


Z)0/?*0, DO - 
DOi?*0, DO . 


K-TT+ 




11 7+'*-'' 

13 4:+^-'^ 


3.47 
3.34 


5.0li;|±0.6 
5.3li;8 ±0.6 


4.3o- 
3.60- 


BO ^ 


£)0K*0^ DO - 


K-tt+ttO 


Q q+6.4 


16.715-5 


3.34 


3.512;? 


1.7(7 


BO ^ 


£)0k*o, simultaneous fit 




41 -oil.! 


10.15 


4.8li;J ±0.5 


5.6cr 


bO ^ 


Z)*0_K'Oj simultaneous fit 




r, y + 3.0 

^•'-2.4 


1.98 


< 6.6 90% CL 


1.4cr 


BO ^ 


D*0K*0^ simultaneous fit 


"■^-4.5 




2.68 


< 6.9 90% CL 


1.4cr 


bO ^ 


D0k*0 , simultaneous fit 


1-41?;^ 


9.21-^ 


10.15 


< 1.8 90% CL 




BO ^ 


D*OK*o ^ simultaneous fit 


i-2tlJ 


0-oir, 


2.68 


< 4.0 90% CL 
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FIG. 3^ AE distributions for the B° d*°K'-*'>° and 
BO i:)(*)Ox*o candidates. Open histo grams represent the 
experimental data and curves show the results of the fits. 



and we set 90% confidence level (CL) upper limits for 
these final states. We do not observe a significant sig- 
nal for the f)i*)Oj^*o decays and also present up- 
per limits for them. Figure |3I shows the AE distribu- 
tions for B° D*°K^*')" and B° d{*)^K*^ candi- 
dates. The upper limit N is calculated from the relation 
Jp^ C{n)dn = 0.9 C{n)dn, where £(n) is the maxi- 
mum likelihood with the signal yield equal to n. We 
take into account the systematic uncertainties in these 
calculations by reducing the detection efficiency by one 
standard deviation. 

As a check, we apply the same analysis procedure to 
B° D+[K'^K+]Tr- and B° D*+[D°Tr+]K- decay 
chains. The estimated branching fractions of B{B'^ — > 
D+TT-) = (2.5 ± 0.3) X lO-^and B{B° D*+K-) ^ 
(1.7 ±0.2) X 10~'*(statistical errors only) , are consistent 
with the world average values 

The following sources of systematic errors are found to 
be significant: tracking efficiency (2% per track), kaon 
identification efficiency (2%), 7r° efficiency (6%), Kg re- 
construction efficiency (6%), efficiency for slow pions 
from D*° -> D°TT° decays (8%), D'^*^° branching frac- 
tion uncertainties (2% - 6%), signal and background 
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shape parameterization (4%) and MC statistics (2% - 
3%). The tracking efficiency error is estimated using ry 
decays to 77 and tt+tt^tt'^. The kaon identification uncer- 
tainty is determined from £)*+ D°7r+, — > K~Tr+ 
decays. The tt*^ reconstruction uncertainty is obtained 
using D'^ decays to K^tt^ and if^7r+7r°. We assume 
equal production rates for B^B~ and B'^B'^ pairs and 
do not include the uncertainty related to this assump- 
tion in the total systematic error. The overall systematic 
uncertainty is found to be f f% for B° D°K^*'>° and 
14% for B° D*^K''*)°. 

In summary, we report the first observation of B^ — > 
£)Oj^(*)o (Jecays. The branching fractions B{B'^ 
D°K°) = {5.0tli ± 0.6) X 10-5 and 3(13° D"K*") = 
{i.8tl;l±0.5) X 10- are measured with 5.1a and 5.6a 
statistical significance, respectively. Note that we ignore 
the possible contribution of B'^ — > D^K^ to the former 
result, since we do not distinguish between K'^ and . 
No significant signal is observed in the — > D*'^K''*^^ 
final states. The corresponding upper limits at the 90% 
CL a,re B{B° D*°K°) < 6.6 x 10"^ and B{B° 
D*°K*") < 6.9 x 10-5. We also set the 90% CL up- 
per hmits for the Vub suppressed B° £)(*)o^*o decays: 
BiB° D°K*°) < 1.8 X 10-5 and B{B° ^ £)*o^*0) ^ 
4.0 X 10-5. 
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